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ABSTRACT 

Context. X-ray observations unveiled the existence of enigmatic point-like sources at the centre of young supernova remnants 
(SNRs). These sources, dubbed Central Compact Objects (CCOs), are thought to be neutron stars formed by the supernova 
explosion. However, their multi-wavelength phenomenology is surprisingly different from that of most young neutron stars. 
Aims. The aim of this work is to understand the nature of the CCO 1WGA J1713.4-3949 in the G347.3-0.5 SNR, through deep 



£ — , optical and IR observations, the first ever performed for this source. 
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Methods. By exploiting its derived Chandra X-ray position we carried out optical (BVI) observations with the NTT and Adaptive 
Optics IR (JHK 3 ) observations with the VLT. 

Results. We detected two faint (I« 23.5, I« 24.3.) patchy objects in the NTT images, close to the Chandra error circle. They 
were clearly resolved in our VLT images which unveiled a total of six candidate counterparts (17.8 < H < 20.3) with quite red 
colours (H-Ks~0.6). If they are stars, none of them can be associated with 1WGA J1713.4— 3949 for the most likely values of 
distance and hydrogen column density. The identification of the faintest candidate with the neutron star itself can not be firmly 
excluded, while the identification with a fallback disk is ruled out by its non-detection in the I band. No other candidates are 
^> ! detected down to B~ 26, V~ 26.2, I~ 24.7, H~21.3 and K~20.5. 

■ Conclusions. Our high-resolution IR imaging of unveiled a few objects close/within the Chandra X-ray position of 1WGA 
' J1713.4— 3949. However, at present none of them can be firmly identified as its likely counterpart. 

Key words. Stars: neutron, Stars: individual: 1WGA J1713.4-3949 



1. Introduction 

X-ray observations have unveiled the existence of pecu- 
liar classes of Isolated Neutron Stars (INSs) which stand 
apart from the family of more classical radio pulsars in 
being radio-silent and not powered by the neutron star 
rotation but by still poorly understood emission mecha- 
nisms. Some of the most puzzling classes of radio-silent 
INSs are identified with a group of X-ray sources detected 
at the centre of young (~ 10-40 kyears) supernova rem- 
nants (SNRs), hence dubbed Central Compact Objects or 
CCOs (Pavlov et al. 2002). 

Send offprint requests to: R. P. Mignani 

* Based on observations collected at the European Southern 
Observatory, Paranal, Chile under programme ID 073. D- 
0632(A),077.D-0764(A) 



Although the SNR associations imply ages of the or- 
der of a few kyears, their X-ray properties make CCOs 
completely different from the other young INSs in SNRs 
(Pavlov et al. 2004; De Luca 2008). Oly two of them ex- 
hibit X-ray pulsations, with periods in the ~ 100-400 ms 
range, and the measured upper limits on the period deriva- 
tives yield spin down ages > 10 3 exceeding the SNR age. 
Furthermore, their X-ray spectra are not purely magne- 
tospheric but have strong thermal components. Finally, 
they are not embedded in pulsar wind nebulae (PWN). 
The discovery of long term X-ray flux variations (Gotthclf 
et al. 1999) and of a 6.7 hours periodicity (e.g. De Luca 
et al. 2006) in the RCW 103 CCO further complicated 
the picture, suggesting cither a binary system with a low- 
mass companion, or a long-period magnetar (De Luca et 
al. 2006; Pizzolato et al. 2008). For other CCOs, the in- 
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voked scenarios involve low-magnetised INSs surrounded 
by debris disks formed after the supernova event (Gotthelf 
& Halpern 2007; Halpern et al. 2007), isolated accreting 
black holes (Pavlov et al. 2000), and dormant magnetars 
(Krause et al. 2005). In the optical/IR, deep observations 
have been performed only for a handful of objects (see De 
Luca 2008 for a summary) but no counterpart has been 
identified yet, with the possible exception of the Vela Jr. 
CCO (Mignani et al. 2007a). 

One of the CCOs which still lack a deep opti- 
cal/IR investigation is 1WGA J1713. 4-3949 in the young 
(< 40 kyears) G347.3-0.5 SNR. Discovered by ROSAT 
(Pfeffermann & Aschenbach 1996), the source was re- 
observed with ASCA and identified as an INS due to 
its high-temperature spectrum and the lack of an optical 
counterpart (Slane et al. 1999). 1WGA J1713.4-3949 was 
later observed with RXTE, Chandra and XMM (Lazendic 
et al. 2003; Cassam-Chenai' et al. 2004), with all the ob- 
servation consistent with a steady X-ray emission. The 
X-ray luminosity is L .5-iofceV ~ 6 1 34 (d/6 kpc) 2 erg 
s^ 1 , where 6 kpc is the originally estimated SNR distance 
(Slane et al. 1999). A revised distance of 1.3 ±0.4 kpc was 
recently obtained by Cassam-Chenai et al. (2004) . The X- 
ray spectrum can be fitted either by a blackbody, likely 
produced from hot polar caps, plus a power-law (kT ~ 
0.4 keV; r - 4; N H ~ 9 10 21 cm" 2 ), or by two black- 
bodies (fcTi - 0.5 keV; kT 2 ~ 0.3 keV; Njj ~ 5 10 21 
cm -2 ). No X-ray pulsations were detected so far (Slane 
et al. 1999; Lazendic et al. 2003), nor any radio counter- 
part (Lazendic et al. 2004) , thus strengthening the case for 
1WGA J1713.4-3949 to be a member of the CCO class. 

Here we present the results of the first, deep optical/IR 
observations of 1WGA J1713. 4-3949 performed with the 
ESO telescopes. Observations are described in Sect. 2, 
while the results are described and discussed in Sect 3 
and 4, respectively. 



2. Observations and data reduction 

2.1. Optical observations 

1WGA J1713.4-3949 was observed on June 13th 2004 at 
the ESO La Silla Observatory with the New Technology 
Telescope (NTT). The telescope was equipped with the 
second generation of the SUperb Seeing Imager (SUSI2). 
The camera is a mosaic of two 2000x4000 pixels EEV 
CCDs with a 2x2 binned pixel scale of 0'.'16 (5'. 5 x 5'. 5 
field of view). 

Repeated exposures were obtained in the broad-band 
B, V, and I filters. The SUSI2 observations log is sum- 
marised in the first half of Table [TJ Observations were 
performed with the target close to the zenith and under 
reasonably good seeing conditions (~ 1"). Since the tar- 
get was always centred on the left chip, no dithering was 
applied to the B and V-band exposures while the I-band 
ones were dithered to compensate for the fringing pat- 
tern affecting the CCD at longer wavelengths. Both night 
(twilight flat fields) and day time calibration frames (bias, 



dome flat fields) were acquired. Unfortunately, due to the 
presence of clouds both at the beginning and at the end of 
the night no standard star observations were acquired. As 
a reference for the photometric calibration we then used 
the closest in time zero points regularly measured using 
Landolt stars (Landolt 1983) as part of the instrument cal- 
ibration plan and available in the photometry calibration 
database maintained by the NTT/SUSI2 team. According 
to the zero point trending plots □, we estimate a conser- 
vative uncertainty of ~ 0.1 magnitudes on the values ex- 
trapolated to the night of our observations. 



2.2. Infrared observations 

1WGA J1713.4-3949 was observed on May 23rd and 24th 
2006 at the ESO Paranal Observatory with NAos COnica 
(NACO), the adaptive optics (AO) imager and spectrom- 
eter mounted at the VLT Yepun telescope. In order to 
provide the best combination between angular resolution 
and sensitivity, we used the S27 camera with a pixel scale 
of 0'.'027 (28" x 28" field of view). As a reference for the 
AO correction we used the GSC-2 star S230012111058 
(V = 14.3), positioned 11"5 away from our target, with 
the VIS (4500 — 10000 A) dichroic element and wavefront 



Observations were performed in the H and K s bands. 
To allow for subtraction of the variable IR sky background, 
each integration was split in sequences of short randomly 
dithered exposures with Detector Integration Times (DIT) 
of 24 s and 5 exposures (NDIT) along each node of the 
dithering pattern. The NACO observations log is sum- 
marised in the second half of Table [TJ For all observa- 
tions, the seeing conditions were on average below ~ 0'.'8. 
Unfortunately, since the target was always observed at the 
end of the night, the airmass was always above 1.4. Sky 
conditions were photometric in both nights. On the first 
night, the second and third K s -band exposure sequence 
were aborted because the very high airmass prevented 
Because of their worse image quality and their much lower 
signal-to-noise, these data are not considered in the fol- 
lowing analysis. The K s -band exposure sequence obtained 
on the second night was interrupted despite of the very 
good seeing because of the incoming twilight. Thanks to 
the combination of good seeing and low airmass, the H- 
band exposure is the one with the best image quality. 
Night (twilight flat fields) and day time calibration frames 
(darks, lamp flat fields) were taken daily as part of the 
NACO calibration plan. Standard stars from the Persson 
et al. (1998) fields were observed in both nights for pho- 
tometric calibration. 



http:/ /www.ls. eso.org/lasilla/sciops/ntt/susi/docs/susiCounts. html 
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Table 1. Log of the NTT/SUSI2 (first half) and 
VLT/NACO (second half) observations of the 1WGA 
J1713.4— 3949 held. Columns report the observing epoch, 
the filter, the total integration time, and the average see- 
ing and airmass. 



yyyy-mm-dd 


Filter 


T(b) 


Seeing (") 


Airmass 


2004-06-13 


B 


3200 


1.14 


1.07 


2004-06-13 


V 


6400 


1.12 


1.03 


2004-06-13 


I 


3150 


1.0 


1.16 


2006-05-24 


Ks 


1800 


0.66 


1.45 


2006-05-24 


Ks 


360 


0.95 


1.70 


2006-05-24 


Ks 


600 


0.78 


1.81 


2006-05-25 


H 


2400 


0.62 


1.42 


2006-05-25 


Ks 


1200 


0.40 


1.74 



2.3. Data reduction 

The NTT/ SUSI2 data were reduced using standard rou- 
tine available in the MIDAS data reduction package^. After 
the basic reduction steps (hot pixels masking, removal 
of bad CCD column, bias subtraction, flat held correc- 
tion), single science frames were combined to filter cos- 
mic ray hits and to remove the fringing patterns in the 
I-band. The astrometry was computed using the coordi- 
nates and positions of 61 stars selected from the 2MASS 
catalogue (Skrutskie et al. 2006). For a better compari- 
son with the VLT/NACO IR images, the I-band image 
was taken as a reference. The pixel coordinates of the 
2MASS stars (all non saturated and evenly distributed 
in the field) were measured by fitting their intensity pro- 
files with a Gaussian function using the GAIA (Graphical 
Astronomy and Image Analysis) tooH. The ht to celes- 
tial coordinates was computed using the Star link pack- 
age ASTR0I"Q. The rms of the astrometric fit residuals was 
w 0'.'09 per coordinate. After accounting for the 0'.'2 con- 
servative astrometric accuracy of 2MASS (Skrutskie et 
al. 2006), the overall uncertainty to be attached to our 
astrometry is finally 0'.'24. 

The VXT data were processed through the ESO NACO 
data reduction pipelinqj. For each band, science frames 
were reduced with the produced master dark and flat held 
frames and combined to correct for the exposure dither- 
ing and to produce cosmic-ray free and sky-subtracted 
images. The photometric calibration was applied using 
the zero point provided by the JVACO pipeline, computed 
through fixed aperture photometry. The astrometric cali- 
bration was performed using the same procedure described 
above. However, since only five 2MASS stars are identified 
in the narrow NACO S27 camera field of view, we com- 
puted the astrometric solution using as a reference a set 

2 http://www.eso.org/sci/data- 
processing/ software / esomidas / 

3 star-www.dur.ac.uk/ pdraper/gaia/gaia.html 

4 http:/ /star-www. rl.ac.uk/Software/software. htm 

5 www.eso.org/observing/dfo/quality/NACO/pipeline 



of 23 secondary stars found in common with the SUSI2 
I-band image, calibrated using 2MASS. The rms of the 
astrometric fit residuals was then « 0'.'06 per coordinate. 
By adding in quadrature the rms of the astrometric fit 
residuals of the SUSI2 I-band image and the average as- 
trometric accuracy of 2MASS we thus end up with an 
overall accuracy of 0'.'25 on the NACO image astrometry. 

3. Data analysis and results 

3.1. Astrometry 

We derived the coordinates of 1WGA J1713.4-3949 
through the analysis of unpublished Chandra observa- 
tions. The field of 1WGA J1713.4-3949 was observed 
on April 19th 2005 with the ACIS/I instrument for 9.7 
ks. Calibrated (level 2) data were retrieved from the 
Chandra X-ray Center Archive and were analysed using 
the Chandra Interactive Analysis of Observations 
software (CI AO v3.3). In order to compute the target po- 
sition, we performed a source detection in the 0.5-10 keV 
energy range using the wavdetect task. The source co- 
ordinates turned out to be a(J2000) = 17 ,l 13 m 28.32 s , 
5(72000) = -39°49'53'/34, with a nominal uncertainty of 
- 0"8 (99% confidence levelf|. The identification of a field 
X-ray source with the bright star HD 322941 at a posi- 
tion consistent with the one listed in the Tycho Reference 
Catalog (H0g et al. 2000) confirmed the accuracy of 
the nominal Chandra astrometric solution. Unfortunately, 
since no other field X-ray source could be unambiguously 
identified with catalogued objects, it was not possible to 
perform any boresight correction to the Chandra data in 
order to improve the nominal astrometric accuracy. 

The computed 1WGA J1713.4-3949 position is shown 
in Fig. 1, overplotted on the NTT/SUSI2 I-band and on 
the VLT/NACO H-band images. In the former (Fig.l-left) 
a faint and patchy object is clearly detected northeast of 
the Chandra error circle (1= 23.5 ± 0.3) and a fainter one 
(1= 24.3 ± 0.4) is possibly detected south of it. However, 
in both cases their patchy structure makes it difficult to 
determine whether their are single, or blended with unre- 
solved field objects. No other object is detected within or 
close to the Chandra error circle down to B~26, V^26.2 
and I~24.7 (3ct). However, due to the better seeing condi- 
tions (see Table [1]) and to the sharper angular resolution, 
five objects are clearly detected in the VLT/NACO im- 
age (Fig. 1-right). Of these, object 413 falls within the 
Chandra error circle. A sixth fainter object (479) is possi- 
bly detected, albeit at very low significance. They are all 
point-like and compatible with the on-axis NACO PSF. 
Objects 401 and 403 are identified with the two faint ob- 
jects detected in the NTT/ SUSI2 I-band image northeast 
and southeast of the Chandra error circle, respectively. 
The former might be actually a blend of objects 401 and 
400, whose angular separation (« 0'.'6) is smaller than the 
PSF of the NTT/SUSI2 image. We thus take the measured 

6 http : / /exc .harvard .edu/cal / ASPECT/celmon/ 
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Fig. 1. (left) 20"x20"JVTT/S!7SI2I-band of the 1WGA J1713.4-3949 field, (right) 6"x6"VLT/NACO H-band image 
of the same field. The area corresponds to the square overplotted on the left hand side image. The circle marks the 
Chandra position of the X-ray source, where the radius (0V85) accounts both for the intrinsic absolute accuracy of the 
Chandra coordinates (99% confidence level) and for the uncertainties of our astrometric calibration (see Sect. 2.3). 
Possible counterparts are labelled on the right hand side image, with object 479 detected only at < 5o\ 



Table 2. VLT/NACO H and K s -band photometry 
and colour of the candidate counterparts of 1WGA 
J1713.4-3949. 



ID 


H 


K s 


H-K s 


380 


19.75 ± 0.14 


19.32 ± 0.11 


0.43 ± 0.18 


400 


18.98 ± 0.13 


18.49 ± 0.09 


0.49 ± 0.16 


401 


17.82 ± 0.13 


17.32 ± 0.08 


0.50 ± 0.15 


403 


18.47 ± 0.13 


17.87 ± 0.08 


0.60 ± 0.15 


413 


18.63 ± 0.13 


18.31 ± 0.09 


0.33 ± 0.16 


479 


20.34 ± 0.15 


19.61 ± 0.12 


0.73 ± 0.20 



magnitude (1= 23.5 ± 0.3) of object 401 with caution. All 
the objects detected in the NACO H-band image are also 
detected in the longest 1200 and 1800 s K s -band ones (see 
Table [l}. No other object is detected close to the Chandra 
error circle down to H~21.3 and K s ~20.5 (3c). 

3.2. Photometry 

We computed objects magnitudes in the NACO im- 
ages through PSF photometry using the suite of tools 
Daophot (Stetson 1992) and applying the same procedures 
described in Zaggia et al. (1997) and applied in Mignani 
et al. (2007a) and De Luca et al. (2008). Since the NACO 
PSF is largely oversampled, we re-sampled the images 
with a 3 x 3 pixels window using the swarp prograrrQ 
to increase the signal-to-noise ratio. As a reference for 

7 http://terapix.iap.fr/ 



our photometry we used the co-added and re-sampled ri- 
band image to create a master list of objects, which we 
registered on the K s -band one and used as a mask for 
the object detection. For each image, the model PSF was 
calculated by fitting the profile of a number of bright but 
non-saturated reference stars in the field and used to mea- 
sure the objects fluxes at the reference positions. Our pho- 
tometry was calibrated using the zero points provided by 
the NACO pipeline after applying the aperture correction, 
with an attached errors of ~ 0.13 and ~ 0.08 magnitudes 
in H and K s , respectively. Single band catalogues were 
then matched and used as a reference for our colour analy- 
sis. The IR magnitudes of our candidates are listed in Tab. 
[2j We used the K s -band photometry performed on the two 
consecutive nights to search for variability on time scales 
of hours. However, none of our candidates shows flux varia- 
tions larger than 0.1 magnitudes, which is consistent with 
our photometric errors. Fig. 2 shows the H,H-K S colour 
magnitude diagram (CMD) for our candidates as well as 
for all objects detected in the field. None of the candi- 
dates is characterised by peculiar colours with respect to 
the main sequence of the field stellar population, which 
suggest that they are main-sequence stars. 



4. Discussion 

To determine whether one of the detected objects is the 
IR counterpart to 1WGA J1713. 4-3949, we investigated 
how their observed properties fit with different scenarios. 
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Fig. 2. H,H-K S CMD of all stars detected in the 
VLT/NACO field. All candidates identified in Fig. 1 are 
marked in red and labelled accordingly. No interstellar ex- 
tinction correction has been applied. 

4.1. A binary system 

If our candidates are stars, we considered the possibil- 
ity that one of them is the companion of the 1WGA 
J1713.4— 3949 neutron star. Their observed colours are 
quite red (0.4< H-K s <0.7), which suggests that they 
might be intrinsically red late- type stars. To be compati- 
ble with the observed range of H-K s (Ducati et al. 2001), 
e.g. an M-type main sequence star should be reddened by 
an amount of interstellar extinction corresponding to an 
N H ~ 10 22 cm~ 2 (Predhel & Schmitt 1995). This value 
is compatible with the largest values obtained from the 
spectral fits to 1WGA J1713.4-3949 (Lazendic et al. 2003; 
Cassam-Chena'i et al. 2004). For the originally proposed 
1WGA J1713.4-3949 distance of 6 kpc (Slane et al. 1999) 
an M-type star with such an high extinction should be 
at least ~0.7 magnitudes fainter than our faintest can- 
didate (object 479). An early to mid M-type star would 
be compatible with the revised distance of 1.3 ± 0.4 kpc 
(Cassam-Chena'i et al. 2004) but it would be detected in 
our NTT/SUSI2 image at I- 20.2 - 21.7. Thus, we con- 
clude that if our candidates are stars none of them can 
be associated with 1WGA J1713.4-3949. Our optical/PR 
magnitude upper limits only allow an undetected compan- 
ion of spectral type later than M. 

4.2. An isolated neutron star 

If 1WGA J1713.4-3949 is indeed an INS, we can then 
speculate if one of our candidates is the neutron star it- 
self. Due to the paucity of the neutron stars observed in 
the IR (e.g. Mignani et al. 2007b) and to the lack of well- 
defined spectral templates, it is very difficult to estimate 



their expected IR brightness. This is even more difficult 
for the CCO neutron stars, none of which has been un- 
ambiguosuly identified so far (e.g. De Luca 2008). In the 
best characterised case of rotation-powered neutron stars 
one can deduce that the magnetospheric IR and X-ray lu- 
minosities correlate (Mignani et al. 2007b; Possenti et al. 
2002). By assuming, e.g. a blackbody plus power law X-ray 
spectrum for 1WGA J1713.4-3949 (Lazendic et al. 2003; 
Cassam-Chena'i et al. 2004) we then scaled the magneto- 
spheric IR-to-X-ray luminosity ratio of the Vela pulsar, 
taken as a reference because of its comparable age (~ 10 
kyears). After accounting for the corresponding interstel- 
lar extinction we thus estimated K s ~ 19.7 for 1WGA 
J1713.4— 3949, i.e. similar to the magnitude of object 479 
(K s ~ 19.6). Since also the magnetospheric optical and 
X-ray luminosities correlate (e.g., Zharikov et al. 2004), 
we similarly estimated Bw 28.3 for 1WGA J1713.4-3949 
which, however, is below our NTT/ SUSI2 upper limit 
(B>26). Thus, a neutron star identification can not be 
firmly excluded. 

4.3. A fossil disk 

As discussed in Sect.l, some CCO models invoke low- 
magnetised INSs surrounded by fallback disks. So, the 
last possibility is that we detected the IR emission from 
such a disk. We note that the IR-to-X-ray flux ratio for 
1WGA J1713.4-3949 would be w lO" 3 - 10" 2 , i.e. much 
larger than that estimated for the anomalous X-ray pul- 
sar 4U 0142+61 (Wang et al. 2006), the only INS with 
evidence of a fallback disk. However, we can not a priori 
rule out the fallback disk scenario. We computed the puta- 
tive disk emission using the model of Perna et al. (2000), 
which accounts for both for the contribution of viscous 
dissipation as well as that due to reprocessing of the neu- 
tron star X-ray luminosity. As a reference, we assumed 
the X-ray luminosity derived for the updated distance of 
1.3 ± 0.4 kpc (Cassam-Chena'i et al. 2004). For a nomi- 
nal disk inclination angle of 60° with respect to the line 
of sight, the unknown model parameters are the disk in- 
ner and outer radii (i?i n ,i? ut), and accretion rate (M). 
We thus iteratively fitted our data for different sets of the 
model parameters. For the dimmest candidate we found 
that the IR fluxes would be consistent with a spectrum 
of a disk (R in = 0.28i? Q , R out — 1.4R©) whose emission 
is dominated by the reprocessed neutron star X-ray lu- 
minosity (Fig. 3), similarly to the case of 4U 0142+61. 
However, such a disk should be detected in the I band, 
with a flux ~ 1.5 magnitude above our measured upper 
limit, as shown in Fig. 3. The overprediction of the optical 
flux is even more dramatic for a disk that fits the brighter 
counterparts. We thus conclude that, if the neutron stra 
has a disk, it was not detected by our observations. 

5. Conclusions 

We performed deep optical and IR observations of the 
CCO 1WGA J1713.4-3949 in the G347.3-0.5 SNR, the 
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Fig. 3. Dereddened IR spectra of the 1WGA 
J1713.4— 3949 candidate counterparts. Dotted lines 
are drawn for guidance. The BVI bands upper limits are 
indicated. The solid line is the best fitting disk spectrum 
(i? in = O.28i? , Ro Ut = IAR G ) for object 479. 

first ever performed for this source, with the JVTT and 
the VLT. We detected a few objects close to the derived 
Chandra X-ray error circle. However, if they are stars the 
association with the CCO would not be compatible with 
its current values of distance and hydrogen column den- 
sity. Similarly to the cases of the CCOs in PKS 1209-51, 
Puppis A (Wang et al. 2007), Cas A (Fesen et al. 2006) and 
RCW 103 (De Luca et al. 2008), our results argue against 
the presence of a companion star, unless it is later than 
M-type, and favour the INS scenario. The identification of 
the faintest candidate with the neutron star itself can not 
be firmly excluded, while the identification with a fallback 
disk is ruled out by its non-detection in the I band. Thus, 
we conclude that the 1WGA J1713. 4-3949 counterpart 
is still unidentified. Deeper optical/IR observations are 
needed to pinpoint new candidates. Although the source 
is apparently steady in X-rays, flux variations as observed 
in the RCW 103 CCO (Gotthelf et al. 1999) can not be 
a priory excluded. A prompt IR follow-up would then in- 
crease the chances to identify the 1WGA J1713. 4-3949 
counterpart. 
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